Bi-partite mode entanglement of bosonic condensates on tunnelling graphs 
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We study a set of L spatial bosonic modes localized on a graph T. The particles are allowed to tunnel from 
vertex to vertex by hopping along the edges of T. We analyze how, in the exact many-body eigenstates of the 
system i.e., Bose-Einstein condensates over single-particle eigenfunctions, the bi-partite quantum entanglement 
of a graph vertex with respect to the rest of the graph depends on the topology of V. 
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The possibility of exploiting the quantum features of 
bosonic particles e.g., cold bosonic atoms, living on coupled 
spatial lattices to the aim of Quantum Information Process- 
ing (QIP) [Jl|] has been recently addressed in the literature 
[H pi H, B]. These systems provide also a unique opportunity 
to investigate fascinating coherent phenomena e.g., quantum- 
phase transitions 

In this note we shall study a simple problem related to this 
more general context. We shall consider a set of N bosonic 
particles hopping between the L vertices of a graph T, we will 
assume the on-vertex self-interaction terms to be zero. The as- 
sociated elementary quadratic Hamiltonian is exactly solvable 
and many-body eigenstates are simply given by Bose-Einstein 
consensates (BECs) over single-particle wavefunctions. This 
kind of abstract situation could be realized, for istance, in a 
optical lattice loaded with cold atomic atoms that can tun- 
nel from different local traps and with atom self-interactions 
somehow switched off [|| . 

The aim is to analyze the role of the graph topology in 
determining, in those many-body eigenstates, the bi-partite 
quantum entanglement of a vertex with respect to the rest of 
the graph vertices. In particular one can address the issue of 
bi-partite entanglement in the ground-state of the system and 
how e.g., for QIP purposes to optimize it by graph designing 
(for a related study see [0 [§ ||). 

It is worthwhile to stress that in this paper the view of quan- 
tum entanglement in system of indistinguishable particles is 
the one based on modes advocated in Refs. [|l0| [HI [j~2|] rather 
than the complementary one based on particles [fl3|] . 

Let us start by recalling the basic kinematical framework an 
to lay down the basic notations. The quantum state-space as- 
sociated with graph T is given by the tensor product of L linear 
oscillator Fock spaces fir — &>j£rsp&n{\nj}}^. =0 . Since we 
are mostly interested in massive particles e.g., atoms, we will 
focus on sectors of Hr with definite total particles number 



\C(m,n 2 ,...,n L )\' 



(2) 



The prime in the above sum simply reminds that the condition 
X^=2 = N — m must be fulfilled. The crucial, though ob- 
vious thing, to notice here is that the constraint of fixed total 
particle-number results in a diagonal reduced density matrix, 
and that such a matrix can be always seen as an operator over 
the finite-dimesional space C with M > N. This remark re- 
lieves us to face with the subtleties of entanglement definition 
in truly infinite-dimensional spaces [ fj~4| ] 
Let the hamiltonian be 



H[A] 
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where 1) the 6i's are bosonic modes, 2) A :— (Aij)ij £ 
Ml(Z 2 ), {Z 2 = {0, 1}) is an symmetric matrix. We will con- 
sider the case in which A is an adiacency matrix of a graph 
T = (V, E), where V = {1, . . . , L} is the set of vertices 
and E is the set of edges, £ E iff Ay ^= 

By diagonalizing A one gets H[A] = X)fc=i ^kB' k Bk, 
where u>k are the A-eigenavalues and B\ = X^=i Ukjbj are 
new bosonic modes (U £ Ml(<C) is unitary). 

Let us now consider a non-degenerate eigenvalue uj\ 
of A and a N particles condensate over it. If B\ = 
(U~n, U12, ■ ■ ■ , U\l) denotes the associated eigenvector, one 
has 
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The reduced density matrix associated to the ith mode is given 
by 



Given a state I*) £ H 
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we are here interested 



to the on-site reduced density matrix; if |^) 
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vertex 
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one has, say for the first 
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-Tr^lBfXBf I = ^ p^(B 1 )\m){m\ (5) 
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where 
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where 

S*(i,m) = {(j„) G IN£7#{j n = i}=m} 

L 

m) := {(nj) G IN^r / m = N, m = m}. (7) 
i=i 

Now, using the fact that X^=i l^ijl 2 = 1 i s not difficult to 
see that one can further rearrange the last expression in Eq.(^) 
in oder to get 

P$(Bi) = (£Q |^ M | 2m (l- |«7i,i|T- m (8) 

This expression is the result we needed. Clearly Eq. (|J) has 
a very simple meaning: the probability p of occupying the 
vertex i (T — {i}) in the single-particle wavefunction B\ is 
given by |C/ H | 2 , (1 - |C/ H | 2 ). Since the BEC over B x is 
the tensor-product of N copies of B\ the probability pfn of 
having m-particle on i is given by a binomial distribution 

f ml P m (^ — p) N ~ m - This classical argument works because 

of the fixed particle-number constraint forces the vertex re- 
duced density matrix to be diagonal i.e., a probability distri- 
bution. 

From now on we will measure entanglement by the von 
Neumann entropy of the reduced density matrix 

e^(B 1 ) := 5(p«( J B 1 )) = -Tr(p«log 2 p«) 

N 

= -$>ml°gaP# (9) 

m=0 

By noting that B\ can be an arbitary single-particle wave- 
function i.e., non necessarily an H[A] eigenstate, one realizes 
that Eq. (^) defines - for any given vertex i of T - a pos- 
itive real-valued function over the single-particle space i.e., 
ejy : C L — * TRq . From expression (||) one readily show that 

• The entanglement of a vertex with respect to the others 
in a BEC depends only on the square amplitude, over 
the considered site, of the single-particle eigenstate we 
are condensing over. 

More formally e§{Wx) = e$(x), (Vz G C L ) for uni- 
taries W belonging to the group U(l) x U(L — 1) (phase 
on the i-th component, arbitary unitary mixing of all the other 
ones). This invariance is, of course, nothing but the invariance 
of entanglement with respect to local transformations. 

• The graph size L does not enter in entanglement proper- 
ties, but possibly through the single particle amplitude 

\Ui,i\. 

(i) 

• It is easy to prove that the functions e N s have a maxi- 
mum for Ui,i = 1 /V2 (see Fig. ^ 




FIG. 1: Entanglement entropy as a function of p := \Un\ 2 for dif- 
ferent particle numbers N. The value has been normalized to the 
maximally available entropy log 2 (7V + 1). 
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FIG. 2: Ratio between max and maximally available entangle- 
ment as a function of the total particle number. 

The highest achieved value for bi-partite mode entangle- 
ment is then given by 

— S-»-r-£(£Wj:) ao) 

m=0 v ' v 7 

this is a monotonic increasing function of N, but - due to 
the well-known properties of binomial coefficients - is mono- 
tonic decreasing fraction of the maximun available entropy 
log 2 (N +1). For N i— ► oo such a fraction seems to attain a 
finite value. We numerically estimated this asimptotic ratio to 
be about 0.57, see Fig (^J). 

Let us define as the set of all undirected graphs with L 
vertices (\Q L \ = 2 i(i ~ 1)/2 ), given r G Ql one has the L 
eigenvectors Bi(T), . . . , Bl(T) of the associated adiacency 
matrix. The real-valued functionals we want to analyse are 
given by 

re0L,— >max fc e$[B fc (r)] (11) 

One could restrict the problem, by considering just the eigen- 
vector B\ (r) associated with the largest eigenvalue of the adi- 
acency matrix of V. This eigenvector corresponds then (see 
Eq. (Q) to the lowest single-particle energy and the BEC |Bf ) 
is a many-body ground state of the Hamiltonian (||). By the 
Perron-Frobenius theorem [ |T3| ] we know that -for connected 
r - B\ is elementwise positive and that the associated eigen- 
value is nondegenerate, hence the ground state is unique. 
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To exemplify this problem let us consider as T the com- 
plete graph minus the diagonal i.e., A[T] = — 
Sij)\i)(j\. By writing this matrix in the following form 

A[r] = L\X)(X\ - 1, \X) := L- 1 / 2 ^^! one imme- 
diately realizes that the A spectrum is given by L — 1 (with 
eigenvector \X)) and by with associated the the L — 1 op- 
erators t>k, (k — l,...,L— 1). The iV-particle ground state 
is therefore provided by putting the N in the k = bosonic 
mode associated with X. The ground-state bi-partite entangle- 
ment is given by (|]) and (||) with \Uij\ = L~ 1/2 . 

If r is a regular graph with connectivity r i.e., all the ver- 
tices have r neighbors, it is fact of elementary spectral graph 
theory Jl5| ] that the highest eigenvalue of the T adiacency ma- 
trix is given by r and the associated eigenvector is given by the 
Fourier mode 1 /yL(l, . . . , 1). Therefore for regular graphs 
maximal bi-partite entanglement is possible just for the dimer 
i.e., L — 2. Notice that for the more general case of one- 
dimensional rings with L (diagonalized by Fourier transfor- 
mation with cyclic boundary conditions) the same value of (||) 
is achieved for all the vertices in all the BECs in single parti- 
cle eigenstates. This fact stems from translational invariance 
which implies that all the single-particle eigenfunctions have 
the same vertex square amplitude i.e., L^ 1 . 

It is interesting to note in passing that the bi-partite graphs 
(V = A U B, (a, b) E E a E A and b E B) the mode 
entanglement associated with BEC over the single-particle 
eigenvalue E is the same as the one associate with eigen- 
value — E. One can realize this fact by performing the fol- 
lowing canonical transformation in the Fock space associated 
with the T modes: Cj — ► (— l) XA ^'Cj, where \A denotes 
the characteristic function of the sub-graph A. One has that 
H[A] — > H[—A]=—H[A], and that the H[A] eigenvectors 
change their components over the Fock basis (g>j-L|n,-) just 
by a phase factor exp(i7r J2jeA n i)- Then the claim follows 
straight away from Eq. (0). Notice also that this symmetry 
property implies that for any initial state \^>) (not necesarily 
an H[A] eigenstate) the on-vertex entanglement dynamics is 
invariant under time-reversal, i.e., S(t) — S(—t), and more- 



over this result holds even in presence of local Hubbard-like 
self interactions Jl6| ] 

For a general number of vertices L, the natural question is? 

What is the graph topology which optimize the on-vertex 
entanglement? 

The answer is not difficult to find out. Let A be the adia- 
cency matrix of the "star" i.e., just the node 1 is connected to 
all the others, Ai j = 8% i. This matrix has two non zero eigen- 
values e± — zt^L — 1 corresponding to the single-particle 
operators 

t±: =>' ± 7I L T^> , ■ (12) 

3=1 

The iV-particle ground state is unique and is given by \b+) = 

0+) tW /VM|0). In view of Eqs (g) and @ the functional 
e]y is maximized for all N by the star graph. Physically this 
means that the star topology optimizes the bi-partite entangle- 
ment in the ground-state BEC. In view of the "monogamy" 
properties of quantum entanglement [|lj this result looks , in a 
sense, rather intuitive. A naive argument is that the star topol- 
ogy is the one with maximall connectivty of the vertex with 
the subgraph with V = (1, . . . , L — 1), this latter in turn is 
totally disconnected and therefore among its vertices there is 
small entanglement. 

In this brief report we studied the mode entanglement in 
Bose-Einstein condensate over a purely tunnel-coupled graph. 
We found an exact expression for such a quantity for arbitray 
graph and particle number. We proved that the star topology 
maximizes the bi-partite entanglement of the spatial mode as- 
sociated to the star center with the rest of the vertices. The 
role of local self-interaction i.e., non-linear terms, as long as 
the practical relevance e.g., implementation, QIP protocols, 
of our abstract though simple analysis is subject of ongoing 
investigations [|6|. 

I thank for valuable inputs, R. Ionicioiu, P. Giorda and R. 
Burioni. 
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